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Introduction: PDC 8.0 includes an implementation 
of a new algorithm to detect and correct step disconti- 
nuities appearing in roughly one of every 20 stellar 
light curves during a given quarter. The majority of 
such discontinuities are believed to result from high- 
energy particles (either cosmic or solar in origin) strik- 
ing the photometer and causing permanent local chang- 
es (typically -0.5%) in quantum efficiency, though a 
partial exponential recovery is often observed [1]. 
Since these features, dubbed sudden pixel sensitivity 
dropouts (SPSDs), are uncorrelated across targets they 
cannot be properly accounted for by the current 
detrending algorithm. PDC detrending is based on the 
assumption that features in flux time series are due 
either to intrinsic stellar phenomena or to systematic 
errors and that systematics will exhibit measurable cor- 
relations across targets. SPSD events violate these as- 
sumptions and their successful removal not only recti- 
fies the flux values of affected targets, but demonstra- 
bly improves the overall performance of PDC 
detrending [1]. 

Algorithm: The algorithm consists of three com- 
ponents: filtering, detection, and correction. The design 
of the LSI detection filter is based on the method of 
Savitzky and Golay [2] with an added multi-scale anal- 
ysis step to improve localization of the peak response. 
Detection and correction are performed on all targets in 
a given channel in an iterative process whereby the 
maximal detector response for each target is evaluated 
and, if an SPSD is detected, a correction is applied. For 
the current iteration, targets in which SPSDs were 
found and corrected are retained in a list and the pro- 
cess is repeated until the list is empty. 

Detection. The detection step begins with precondi- 
tioning of each time series by filling gaps and padding 
endpoints. Conditioned time series are then convolved 
with the detection filter and responses are normalized 
across all targets and cadences on the current channel. 
The maximum normalized detector response for each 
target is evaluated by applying succession of threshold- 
based decisions to determine first whether it is a feasi- 
ble candidate and finally whether it should be labelled 
an SPSD. 

Correction. Correction is a two stage process. The 
first stage estimates a persistent step height from analy- 
sis of the entire flux time series, excluding a short re- 
covery window following the SPSD, which typically 
contains a transient signal. The second stage models 
the recovery window using a series of exponentials of 


varing time constant. The algorithm also detects and 
preserves sinusoids in the signal while removing only 
the step and recovery transient components. 

Results: The algorithm has been tested on (1) Q7 
flight data containing signatures of real SPSD events, 
and (2) Q7 flight data with SPSD-affected targets re- 
moved and 100 simulated SPSDs injected into the time 
series of —12th magnitude targets (11.5 < m < 12.5). 
Detection performance on flight data was estimated at 
P(hit)=0.89 and P(false alarm)=0.007, relative to man- 
ual SPSD identification for a single channel (1912 light 
curves). Performance on simulated SPSDs is shown in 
figure 2 as a function of the magnitudes of injected 
sensitivity drops. 
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Figure 1. Example of a typical discontinuity and its 
correction (channel 13.1, Q7). 
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Figure 2. Detection performance vs. simulated sen- 
sitivity drop (channel 13.1, Q7). 


